BILINEAR EMBEDDING FOR REAL ELLIPTIC DIFFERENTIAL 
OPERATORS IN DIVERGENCE FORM WITH POTENTIALS 



OLIVER DRAGICEVIC AND ALEXANDER VOLBERG 



Abstract. We present a simple Bellman function proof of a bilinear estimate for el- 
liptic operators in divergence form with real coefficients and with nonnegative poten- 
tials. The constants are dimension- free. The p-range of applicability of this estimate 
is (1, oo) for any real accretive (nonsymmetric) matrix A of coefficients. 



1. Motivation 

The bilinear-embedding-type theorems as the one we prove below for general second 
order operators in divergence form with potentials, Lu := —div{AVu)+Vu, where ^ is a 
real accretive matrix and V a nonnegative locally integrable function on R", have so far 
played an important part in estimates of several singular integral operators. In \17\ [TT] 
the bilinear embedding was used for estimates of the Ahlfors-Beurling operator T. 
In particular, that implied a regularity result for solutions of Beltrami equations at 
the critical exponent, which, in its turn, had a geometric significance for the theory of 
quasiregular maps; see [1]. 

When applied to families of operators on spaces of arbitrary dimensions, the bilinear 
embedding tends to give dimension-free estimates with sharp behaviour in terms of p. 
The first such example is [9], where dimension- free estimates were obtained for Riesz 
transforms associated to the classical and the Ornstein-Uhlenbeck Laplacians. This 
approach was continued in [8] in the context of arbitrary Schrodinger operators with 
nonnegative potentials. As a special case, dimension-free estimates for the Hermite- 
Riesz transforms on were confirmed, with linear behaviour with respect to p. 

Recently, A. Carbonaro and the first author [6] studied a class of Riesz transforms 
in the fairly general setting of complete Riemannian manifolds whose Bakry-Emery 
Ricci curvature is bounded from below. The study of such operators was initiated by 
Strichartz |19j . Bakry [5] proved the first dimension- free estimates, whereupon the work 
of Li |12j implied that in the general case the behaviour of the norms with respect to 
p is at most quadratic as p — ?• oo. In j6] this was improved by giving linear estimates. 
They are sharp. Also, while the proofs in [5] and [12j are probabilistic, [6] features the 
first analytic proof. The bilinear embedding was instrumental for obtaining all of these 
results. 

The operators we consider are also called generalized Schrodinger operators, see [7]. 
When A = I we have the usual Schrodinger operators —A -|- V. The fact that A = I 
makes them comparatively easy in the theory of general second order accretive operators 
in divergence form, so the emphasis in, say, [HI [HI El [10] is on independence of estimates 
on the dimension of the underlying manifold. 
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1 



OPERATORS IN DIVERGENCE FORM WITH POTENTIALS 



2 



Yet another place where the biUnear embedding featured (under a shght disguise) is 
the sharp estimate of the weighted Hilbert transform in terms of the Ap norm of the 
weight [TelllH]. 

Of course the dimension-free estimates of Riesz transforms were known before. One 
can find an extensive bibliography in and [12] . We think the approach in [9l [HI [6l [10] 
is a sort of a unified one. Namely, there the proofs are divided into two steps: 1) the 
proof of the bilinear embedding (always dimension-free) and 2) a formula which involves 
holomorphic calculus of the operator L. We would like to emphasize the importance of 
the second step because of its relations with Kato's problem. The holomorphic calculus 
of general second order accretive operators in divergence form also plays crucial part in 
proving estimates for the related Riesz transforms, see e.g. [21 [^. 

Finally, we would like to acknowledge that the completion of this paper was decisively 
motivated by recent appearance of a closely related preprint [3] by Auscher, Hoffman 
and Martell. In particular, [3] explains the connection between the results in this paper 
and square function estimates in the context of accretive matrices. 

Acknowledgements. We are grateful to P. Auscher, S. Hofmann, C. Kenig and J. 
Pipher who advised us on general second order elliptic operators. Special thanks are to 
J. M. Martell for several valuable conversations and remarks. 

We are particularly grateful to Pascal Auscher whose interest made us resume working 
on this paper after a pause of over three years. 



where z ■ w := ZjWj for z,w £ C". 

Let y be a nonnegative locally integrable function on M". We consider the operator 



and the associated semigroup {Pt := e ^^)t^Q. For precise definitions see [3 Section 1.8] 
or [151 Section 4.7]. Denote furthermore f{x,t) = {Ptf){x). 

Next, for a smooth complex-valued function (j) = (j){x,t) on M" x (0,oo) introduce 
\(p\l := |Vi;^p -|- where x = (xi,...,rE„) and V = (dx^, ■ ■ ■ ,dx„) is the spatial 
gradient. That is. 



Take any number p £ [2,oo) and denote by q its conjugate exponent, i.e. 1/p + l/q = 1. 
What follows is the main result of the paper, which we call bilinear embedding theorem. 

Theorem 1. Let the operator L be as above. For any f,g€ C^(M") we have 



2. Elliptic differential operators in divergent form 





n 



10(^,^)1':= E|^(^'*)| +nx)i</>(x,t)r 
j=i 1 



(2) 




\ f{x,t)\^\g{x,t)\^dxdt ^ C^p||/||p|bllg • 
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The theorem is also vahd for 1 < p < 2, in which case the factor on the right becomes 
C-^q. Moreover, it will emerge from the proof of the theorem that one might take 
= Cmax{l,7~^}, where C is some absolute constant. 

When V 
above: 



0, Theorem [T] would obviously follow if we had the following estimate from 



{GLu{x)Ydx ^ C / \u{x)\Pdx 
for all p G (l,oo), where Glu is the square function, i.e. 

Glu{x) := ( r \VPtu{x)\^dt 



(3) 



In general this inequality is proved for p G {q^{L), Even when ^ is a real matrix, 

([3]) holds only for p £ (1,2 + See |2l Corollaries 6.3 - 6.7]. On the other hand, 

an interesting paper by Auscher, Hofmann and Martell [3] shows that the cone square 
function estimate for real matrices A can be made for the whole range of p G (l,oo). 
The proof seems to be rather delicate, and as far as it stands now the estimates depend 
on the dimension. Our theorem has a strange feature of being dimension-free. 

3. Bellman function 

The Bellman function technique was introduced into harmonic analysis by Nazarov, 
Treil and the second author in the 1994 preprint version of their paper |14] . The main 
tool in our proofs will be a particular Bellman function, drawn from a function defined 
by Nazarov and Treil [13]. Their example was later extended by the present authors in 
[8] . In this paper we work with a simplified variant of the function from [8] ; it comprises 
only two variables, while the functions from \13\ [8] have four. 

Throughout this section we assume that p ^ 2, q = p/{p — 1) and 5 = q{q — l)/8 are 
fixed. Observe that 5 ~ {p — 1)~^. 

Let : M+ X M+ — > IR+ be defined by 



uP + 



1 



w 



2 

p \q 

The Bellman function we use is simply the function Q : C x C 



2 2- 
u V 



w 



P>vi . 



Our proofs will to a great extent involve estimates of the first' 
derivatives of Q. Estimating the gradient of (j) is straightforward: 

du^ ^ C{p) max{uP-^,v} and d^4 ^ Cv""^ C{v + 1) . 



(4) 

M_, defined by 

(5) 

and second-order 



(6) 



For the last inequality we used that q — 1 ^ 1. Now the chain rule naturally implies the 
estimates of d(^Q and dr^Q. 

Yet while Q is of class C^, it is not globally C^, because (j) fails to be along 
the curves {v = 0} and {u^ = v'^} in M^. This proved not to be a major obstacle, 
since smoothening Q by standard mollifiers enables one in various contexts to apply Q 
completely rigorously, see [Hill [9]. Hence, for the sake of the clarity of the presentation, 
we will further proceed as if Q were everywhere, referring the reader to [HI |6l [9] for 
the mollifying procedure. A similar smoothening argument was also carried out in |17j . 
With this understanding, the properties of Q are summarized in the theorem below. 
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Theorem 2. For any u,v ^ 0, 

(i) ^ (l){u,v) ^ {l + 6){uP + v^). 

// ^ = (C, ??) G C X C then there exists r = t{\C\, \r]\) > such that 

(ii) -d^Q(0^6{T\dC\^+T-^dv\^) 

(iii) Q(e)-dQ(Oe^5(r|CP+r-i|7?|2). 

The first property follows immediately from the definition of (p, while the last two 
were proved in 0. 

Let us clarify the notation appearing in the preceding theorem. Denote by Zj the 
variables of Q, i.e. 

{zi,Z2) ■■= {C,V) ■ 

Write Zj = tCj + it)j for j = 1,2, so that 

d 1/9 d \ d 1/5 d 

Take a = (cri,cj2) G C^. Then we define 

j=l J J 
As for the second differential form, when ^, cr, ^ € C^, define 

2 / q'^q d^Q d^Q d^Q 

^ \dZkOZrn OZkOZrn OZkOZm OZkOZm 

k,m=l 

Note that this expression is always real. It is also symmetric, in the sense that 

(d2Q(0cT,?> = (d2Q(eX,a> . (8) 
The meaning of ^ above is the inequality 

{-d^Q{Oa,a) ^ 5{T\ai\^ + T~'\a2\^) , 
while (jml) should be understood through d?]), with o" = ^, of course. 
3.1. Setup of the proof. Given functions f,g : M" — C, define 

vi{x,t) = Ptf{x) and V2{x,t) = Ptg{x) (9) 
and V := {vi,V2), i.e. 

vix,t) = {Ptf{x),Ptg{x)). 

Put b := Q o V, that is, 

b{x,t) = Q{Ptf{x),Ptg{x)). (10) 

When proving Theorem [1] we will try to imitate [9]. Let us explain what this means. 

Suppose Tp G C^(]R") is a radial function, V = 1 in the unit ball, -0 = outside the 
ball of radius 2, and < ip < 1 everywhere else. For R > define 

ipnix) := ^(^^ 

Furthermore, let 
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Take R,T > and consider the integral 

lR,Tif,g)-= / ipRix)L'b{x,t)dxdt, 
Jo JK" 

where /, g are the same functions as in (jlU|) . The proof of Theorem [T] will rest on the 
following two estimates of Ir^t (lower and upper, respectively). 

Proposition 1. Let j,5 be as in ([1]) and Q. For every {x,t) G R"' x (0, oo) we have 

L'b{x,t) ^ 26 mm{l, -f}\f {x,t)U\g{x,t)U . 
Proposition 2. For arbitrary test functions f,g and T > 0, 

limsup r [ ijR{x)L'b{x,t)dxdt ^ 11/11^ + ||5||^ 
R^oo Jo JR" 

Proof of Theorem [IJ We use the standard trick of "polarization" . Together the Fatou 
lemma and Propositions [1] and [2] imply 

|/>,OUIg(^,OUrf^^t^ (11/11^ + l|g||^)- 

Now replace f,g by A/, X~^g, respectively, and take the infimum over A > 0. □ 

So what remains is to prove Propositions [1] and [21 This will be done in Sections S] 
andO respectively. 

4. Proof of Proposition [T] 

We would like to apply properties ^ and ([nH) of the function Q that were given in 
Theorem [2j To this end we must first link them to L'b. 

Lemma 1. Fix {x,t) G R^^^ . If [oij] is a real accretive n x n matrix function, then 
L'b{x, t)=f2 ("'^^^('')^(^' + ^(^)[^(^) - dQ{v)v] . (11) 

i,j=l ^ * J / 

Proof. By applying the chain rule repeatedly we calculate 



and 

d'^b I jir\i \ \ ( t ^^^fc ^'3/ d'^Vk 



ffQ(v)zlL HI. \ + ( ^ (v) ^ + —(v) 



dxidxj \ dxi dxj I ^^K^'^k dxidxj dz^ dxidxj 

Now we are ready to compute 



Bb 

L'b = — - div(ylV6) + Vb . (12) 
Fix k and consider only those terms in ()12p which contain dz^Q{v) as a factor. We get 



|£(^^(x,t)) - div(AVi;fc)(x,t) 



which is the same as 



||(t;(x, t)) (L'^;fc)(x, t) - 0) ^(^)^^fc(^, • 
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But each is, by ([9]), a function of the form e~^^ip, therefore L'v^ = 0. 
If we now factor out dz,^Q{v), we get 



where A is the matrix with entries {oij}. Since A was assumed to be real, the term in 
parentheses is again equal to L'vk, hence it disappears. 

Putting together the remaining terms gives (jlip . □ 

Proof of Proposition {1[ We would like to rewrite 

aij{x) ((fQ{v) — {x,t), (13) 

i,j=l * 3 1 

in terms of quadratic forms, i.e. in a way where instead of dx^v and dx -v we would get 
identical vectors. This would clear the road for applying Theorem [2] . 
For any square matrix A = [uij] let 



(14) 



and A = [dij]. If A is real accretive, then so is A, but it is, in addition, also symmetric. 
First note that the symmetry ([8]) quickly implies 

dv , , dv 

a.ij{x) (d^Q{v) 

Consequently, we get 



m=Yl aij{x)U^Q{v) — {x,t),—{x,t) 



k,m=l *J=1 k,m=l *J=1 

+ E si^i-) E «^^(-) 9^ + E ai^(-) E «^.(-) ^ • 

K,m=l «J = 1 k,m=l *J=1 



This is the same as 

k,m=l k,m=l 



Matrix A is real and symmetric. The accretivity condition implies that it is also 
positive-definite, hence it admits a real (positive-definite) square root A^^'^. By intro- 
ducing new vectors 

Gfc := A'/'Vvk=: {9t---,0'j 
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and using that yl^/^Vw^ = A^/'^Vvk we can rewrite the last sum again as 



E 

k,m=l 



^ dzkdzm 

n 2 

E E 

j=l k,m=l 



(V) {Ok ■ @m) 



dzkd 



kOZ„ 



iv) 



^ ^ dzkdz„ 



^ ^ dzkdz„ 



dzkdz. 



This expression we can write in turn as X]"=i {d'^Q{'^)^jj^j)i where 9j := {9j,6j) € 
for j G {1, . . . , n}. So we proved 



We can finally apply Lemma [T] and properties (jn]), (pli|) from Theorem [5J The result is 

n 

L'b{x,t) = {-d'^Q{vWj,Oj) + V{x)[Q{v) - dQ{v)v] 



n 



^ 25 




= 25V^l^iP + |0ilV^b2p + |e2|2. 

Recah that = iVSy?;^ and use that 2^ • ^ = Re (A^ • 0- Therefore, ^ gives 

^ 25min{l,7}Vy|ui|2 + \Vvx\^^JV\v2\^ + jVvap . 
Finally, by ([2|) and (l9|), this is the same as 25min{l,7}|/(x,t)|*|^(x,t)|^, . 

5. Integration by parts: proof of Proposition [2] 



□ 



In order to prove Proposition [2] we will first integrate by parts, as summarized in the 
following inequality. 

Lemma 2. Let A be a real accretive matrix (not necessarily symmetric) and p ^ 2. 
Then 



lim sup 

R-^co Jo 



T 



ipji{x)L'b{x,t) dx dt ^ lim sup 

i?,--s>oo Jo 



Wb\x) -^{x,t)dxdt . 



Proof. We first recall (|12|) and notice that Vb ^ 0. Hence it suffices to prove that 

tpR{x)div{AVb){x,t)dxdt = 0. 



lim sup 

R^oo Jo 



Because of the integration by parts and the uniform boundedness of A it is sufficient to 
show 



lim 

_R— 5>oo 



T 



\V'iPr{x)\ \Vb{x,t)\dxdt = 0. 



Jujn 
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Here ojr := supp Vi/'ij = {x G M" ; ii ^ |x| ^ 2R}. When R is large, 

d(a;_R,supp/) ~ R. 
Since |(V^fl)(x)| ^ i^-lVV'll oo, our inequality boils down to showing 



(15) 



lim 



\Vb{x,t)\ dxdt = 0. 



'0 Juju 

By pop , this will in turn follow once we demonstrate 



lim 



Jujji 



dQ 
dC 



yPtf{x)\ + 



dr] 



VPtg{x)\ dxdt = 0. 



(16) 



For that purpose we resort to the so-called off-diagonal estimates for operators L which 
were proven in [21 Section 2.3]. That is a result stating that for any h, supported in a 
closed set E C M", the following holds: 

||ri/i||i.(P)^Ce-f'^(^'^)'||/i||i2(^). (17) 

Here F C M" is closed and Tt can be anything from Pt,tLPt,\/t\7Pt. While the proof 
in [2] was done for the case V = 0, it can be repeated so as to cover the case of general 
nonnegative potential V G Lj^^{M."'). 

In order to estimate the partial derivatives of Q with respect to ( and i] we first recall 
([6]), which suggests considering separately the domains 

Ai := {(C, r?) G C2 ; \C\P ^ m and A2 := {(C, r?) G ; \C\P > ^ . 

The upper estimates of VQ in Ai and A2 can now be directly inferred from ^ and ([6]). 
Fix t ^ and define, for j G {1, 2}, 

= Vj{t) := {x G M" ; v{x, t) G A^} . 

1.) Let us estimate \d,^Q{v)\\VPtf\ in domain Vi, that is, when \Ptf\^ ^ \Ptg\'^. This 
condition is equivalent to \Ptf\^^^ ^ \Ptg\, therefore, by ([5]) and ([6]), 



dC 



^C\Ptg\. 



(18) 



This and the Holder's inequality imply 



/ 



^{v)\\VPtf\ ^ ^JPtg\\LHu.,)\\VtVPtf\\LH.^) 



By the off-diagonal estimates p7p and the similarity psp we get 

C{f,g,n) _c^2 

Finally, integrate in t and send R to infinity. Now it is clear that 



lim 



/ XViw(^) ^(^;(x,t)) |VPJ(x)|dx(it = 0. 

Jlur Oi, 



2.) In the domain V2 we similarly get 

Applying the Holder's inequality again gives 
dQ, 



I ^v,|^(^')||vpjKc(£ IPjl'^^-^^Y'wvPtfh^ 
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Since the semigroup Pt is L°°-contractive [151 Proposition 4.32] we can continue as 
Note that \ujfi\ = C R"". Together with ()17p this imphes 



C{f,p,n) /2 

which again disappears after performing hm/j_>.oo Jq dt. 
3.) Finahy, the estimates with respect to r/ can be treated without sphtting the cases 
X G Vi,2(t). Indeed, for any (x,t) G M"+^ we have, by and ©, 

^-^{v)\^Cmg\ + l). 

Thus we got a sum of two terms, of which the first is the same as in p8|) . hence it 
become zero after multiplying with jVP^^I, integrating and taking the limit. As for 
the last remaining term (i.e. the constant 1), the Holder's inequality and (|17|) give 



These estimates combine into (|16p . which finishes the proof. □ 
Proof of Proposition In view of Lemma [2] it is enough to prove 

hmsup r [ i;R{xAx,t)dxdt^\\frp+\\g\\l 



We start by another integration by parts: 



T 



db 

ipR{x) — {x,t)dxdt= I ijR{x)[b{x,T) - b{x,0)]dx . 



Remember that Q was a nonpositive function, therefore b{x,T) ^ 0, by (jlOp . Moreover, 
(jlOp and ([5]) also mean we can continue as 

^11 H\f{x)\,\gix)\)i;R{x)dx. 

By Theorem [2] ([i]) we get 



which tends to + \\g\\q as i? — t- oo, by the dominated convergence theorem. □ 
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